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Abstract 

In this article we consider the possibility of observing nonstandard Higgs bosons 
in the H — > ZA — > Zt + t~ channel. We present three benchmark scenarios in the 
NMSSM where H — > ZA is the dominant decay mode for one of the nonstandard 
Higgs bosons while the lightest CP-even Higgs is Standard Model like. Using the latest 
CMS multilepton analysis based on 7 TeV LHC data, we put limits on the signal 
cross-section, which constrain leptophilic scenarios. Projecting to future LHC analyses 
with improvements in background modeling, we show that with O(30) fb _1 of data, 
such a multilepton analysis is very close to constraining our NMSSM benchmarks. As 
we illustrate with a toy model, for light A masses, the large boost of the A makes it 
inefficient to select two hadronic taus, since isolation and the transverse momenta are 
in tension. This efficiency could be improved by including boosted di-tau jets as an 
object in future multilepton analyses. We also discuss different methods to confirm 
this scenario by reconstruction of the run and masses. In particular we consider 
the transverse mass distribution, collinear mass distribution and an analytical solution 
using trial masses. 



1 Introduction 



Recently, ATLAS and CMS announced the discovery of a particle consistent with the Stan- 
dard Model Higgs [U [2] . If the properties of this particle are confirmed to be that of a Higgs 
boson, this will be a major revolution in particle physics, as it would be the first fundamental 
scalar field proven to exist. This discovery opens up the possibility for other fundamental 
scalars and thus motivates searches at the LHC for additional spin zero particles. Moreover, 
in many theories beyond the Standard Model, e.g. supersymmetry or two Higgs doublet 
models, such additional Higgs bosons are required and thus their discovery would be an 
important first step towards uncovering this new physics. 

The phenomenology of such scalars has so far been incompletely explored. In the minimal 
supersymmetric Standard Model, there is a well explored phenomenology of the heavy Higgs 
states H, A with ongoing searches in t + t~ [3j H] and potentially observable bb decays [5]. 
However, scalars can also have extremely challenging signals, many of which were explored 
in the nonstandard Higgs scenario [6]. Thus, it is fruitful to continue to explore promising 
signals for new bosons. 

In this note, we point out an interesting signal with the potential to discover two new 
scalars. The signal process is H — > Z, A — >• £ + £~ ,t + t~ . Such a signal is motivated by 
the recent CMS multilepton search [7j in particular its channels with an onshell leptonic 
Z. Recently, the multilepton search has also been shown to be sensitive to Standard Model 
Higgs modes [8], flavor violating top decays into a Higgs (t — > ch) [9], and two Higgs doublet 
models [10J. These previous analyses did not consider reconstructing hadronic taus and 
focused on an inclusive search using all the CMS multilepton channels. Considering our 
signal's exclusive channels into hadronic taus provides an important additional probe which 
is sensitive to a substantial branching fraction of the signal and also allows a variety of mass 
reconstruction techniques given the limited number of neutrinos. Admittedly, the substantial 
excess in the original CMS multilepton search [TT] in these channels went down in [7j, but 
this signal still remains promising in future updates. 

In this paper, we analyze the prospects a CMS-like multilepton analysis has in discovering 
such a signal. The outline of the rest of the paper is as follows. In Sec. |2j we explore simple 
benchmarks to realize such a signal in the Next-to-Minimal Supersymmetric Standard Model. 
This provides an important existence proof in a motivated theory. More optimistic scenarios 
exist as well in leptophilic two Higgs doublet models. Those interested in the phenomenology 
can skip to Sec. [3j where we estimate the efficiency the CMS multilepton search has on such 
a signal and use them to derive model-independent bounds on the signal rate. In Sec. [4], we 
compare the utility of a variety of mass reconstruction techniques for determining the H, A 
masses. In Sec. |5j we conclude and look at future directions. In Appendix |Aj we give details 
for solving the neutrino momenta given trial masses for the two bosons. 
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2 A simple model for H — > ZA signals 



In this section we discuss the possibility of enhancing the coupling of a CP-even Higgs boson 
(H) to the Z gauge boson and a pseudo-scalar (A). We would like to develop a model 
where the H decays predominantly to ZA and the A in turn decays into t + t~. This model 
has the possibility of explaining the slight excesses observed by the CMS collaboration [7] 
in the 4£ + Or^, 3£ + lr^ and 21 + 2t^ channels with a reconstructed leptonic Z. In these 
channels, the r^'s indicate reconstructed one-prong hadronic decays. Therefore, to realize 
this scenario in a model we need the A to be reasonably light and the branching ratios of 
H — > ZA and A — > t + t~ to be significant. In particular, we are interested in benchmarks 
where m# ~ 200 GeV and ttia ~ 10 — 100 GeV, ensuring that the Z decay is open and the 
H can be produced with reasonable rates. Such a decay is particularly difficult to realize in 
the Minimal Supersymmetric Standard Model due to the decoupling limit constraining both 
the mass of A and the decay of H — > ZA and thus, we have to turn to other theories. 

To demonstrate a realization of such a model we consider the Next-to-Minimal Super- 
symmetric Standard Model (NMSSM) [12]. The superpotential has the form 

W = W Ynk + XH u H d S+^S 3 (1) 

where Wyuk are the usual Yukawa interactions and the hatted fields denote the chiral super- 
fields. The corresponding soft supersymmetry breaking terms are 

VU = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 \S\ 2 + V2 (m x SH u H d - ^S 3 ) . (2) 



Here, we will follow the discussion and notation of Ref. [13J . The relationships to the standard 
NMSSM notation of Ref. [H] are m K = -kA k /\/2 and m A = \A x /\/2. We will work in the 
CP-even Higgs basis (h° v , h° s ) and CP-odd basis A° s ), where 

H° d = ^[(v + K-iG )ce-(H v -iA v )sp] (3) 

Hi = -L[( v + h v + iG )sp-(H v + iA° v )cp] (4) 

S = -L( s + h s + iA% (5) 

v ~ 246 GeV is the Higgs vacuum expectation value (VEV), tan/3 = (v u /v d ), G° is the 
goldstone mode, s is the 5* VEV. An effective /i parameter can be defined /z e fr = \s/V2. As 
h® is rotated in the same way as v, it is the linear combination that gives mass to the W 
and Z and hence is the one that has trilinear couplings to these gauge bosons. Therefore 
the Standard Model-like Higgs boson is the one which has the largest component in the h® 
direction, while the nonstandard Higgs boson is the one that is mostly in the H® direction. 

The minimization conditions can be used to eliminate some of the soft SUSY breaking 
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After substituting these solutions into the scalar potential, there are six remaining free 
parameters A, k, (3,m\,m K and s. In terms of them, the tree-level CP-even mass matrix in 



the (h° v ,H°,h° s ) basis is 
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where 



A 2 
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M 2 



sin 2 2/3, 



i? = -(2Aks + m A ) sin 2/3. 

Similarly, the CP-odd states (A®, A®) have a tree-level mass matrix 

P\2 _ f i^ Ks2 + 2m x s)s^ -v(Xks - m x ) 

A \ -v(\K,s-m x ) (Xk+ ^)v 2 s 2 p + 3sm K 

where the rotation angle 9a satisfies 

4vs(Xks — m\) sin 2/3 



tan 29 a 



v 2 sin 2 2/3(2Aks + mA) - 2s 2 (Xks + 2m x - 3m K sin 2/3) 
In this basis the ZA couplings to the CP-even states have the form 

D^HlD^H u + D»H\D^H d 
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where a& = ad J, - bd u a and in the last line we have rotated into the CP-odd Hi 
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Within the NMSSM, a light pseudo-scalar can exist either in the Peccei-Quinn (PQ) or 
the U(1)r limit. If we take the U (1) ^ axion limit discussed in Ref. [13] . where O(10~ 3 ) < mA, K /f < 
1, the mass of the lightest CP-odd scalar is 
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(GeV) 


(GeV) 


(TeV) 


(GeV) 


BM1 


0.74 


1.10 


1.5 


-9.1 


-10.0 


0.0 


166 


0.5 


BM2 


0.74 


1.10 


1.5 


-9.1 


-10.0 


0.0 


175 


0.5 


BM3 


0.70 


0.78 


1.5 


-4.2 


-40.6 


0.0 


170 


0.5 



Table 1: Model parameters of benchmark scenarios for enhanced H\ — > ZA\ signals. 
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(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 






BM1 


125.9 


263 


470 


9.83 


351 


259 


0.982 


-0.691 


BM2 


125.9 


280 


495 


23.9 


364 


276 


0.985 


-0.685 


BM3 


124.9 


235 


361 


105 


317 


232 


0.995 


-0.657 



Table 2: Higgs mass spectra and normalized coupling to top quark in the benchmark sce- 
narios. 
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final state 


77 


WW* 


ZZ* 


A\A\ 


ZA\ 


A\A\ 


TT 




BM1 


2.4 x lO" 3 


0.18 


2.0 x lO" 2 


8.7 x 10~ 2 


0.77 


0.19 


0.68 


.25 


BM2 


2.6 x 10- 3 


0.20 


2.3 x 10- 2 


1.1 x lO" 4 


0.78 


0.17 


6.4 x 10~ 2 


.023 


BM3 


2.7 x lO" 3 


0.19 


2.1 x 10- 2 


0.0 


0.68 


0.29 


9.0 x 10~ 2 


.026 



Table 3: Relevant branching ratios of the Higgs bosons in the benchmark scenarios and the 
signal rate \x = u (vv^h^za-^Zt+t ) norma ]i zec i ^ Q ^he Standard Model Higgs cross section. 
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tan0„ ~ — - + . (It 



and the rotation angle is 

tan a ~ 

f sin 2/3 

In order to have a light with a large gn°zA\ coupling we require large mixing in the 
CP-odd Higgs sector, 8 a — f • Also, to satisfy charged Higgs limits from top decays [EJUS], 
we have to raise the A2 mass since it is correlated with the charged Higgs. In the R-axion 
limit, the magnitude of m A o is set by Xk{v 2 S2I3 + s 2 ^ 1 ); so the mass constraint implies that 
A and k are both order one. To suppress the g b i H o coupling relative to Qh % za\ we need 
to be at low tan /3. Additionally, maximal mixing in the pseudo-scalar sector along with 



Eq. (18) suggests that s ~ v sin 2/3 ~ v. Large gn°zA° couplings also typically lead to an 
enhancement of the g^o^o^o which reduces the Hi's branching ratio into Standard Model 
decays, which is constrained by the present Higgs boson measurements at the LHC [TJ [2] . 
Hence, i when < m^/2, viable benchmark points need a tuning so that the g^o^o^o 
coupling is suppressed, such that the branching ratios of Hi remain SM-like. 

We note in passing that the PQ-axion is not a useful limit for our benchmarks. In the 



PQ limit with K,m K small, we can neglect the k and m K terms in Eq. (12) and find that 
9a — vr/4 implies s ~ vs2p/2. Since the VEV s ~ v/2 ~ 125 GeV is smaller than in the PQ 
limit, the results small // e // often leads to chargino masses in violation of LEP2 bounds. 

Therefore we consider benchmark points in which the pseudo-scalar is an R-axion and 
use NMSSMtools 3.2.0 [H] to find the benchmark scenarios with NMSSM parameters shown 
in Tab. [T] and physical Higgs boson masses in Tab. [2] As can be seen in the last column of 
Tab. [2J the H\(H2) has production cross sections through gluon fusion ~ 1(1/2) times the 
cross section for a Standard Model Higgs at that mass. The branching ratios are shown in 
Tab. [3} To suppress couplings to down-type fermions we have chosen to only consider 
tan/3 = 1.5 scenarios. The values of A and k are large in order to generate a significant mass 
splitting between the two pseudo-scalar states. Some of these couplings are large enough to 
develop a Landau-pole before the GUT-scale. Hence the UV-completion of such scenarios 
may require Fat-Higgs like models discussed in Ref. [TTj . For these regions of parameter 
space, the Standard Model-like Higgs boson is the H\ state and for A % 1, its tree-level 
mass m^o e ~ 100 GeV. Hence a small amount of stop radiative contributions is needed to 
raise the physical SM-like Higgs to the observed Higgs boson mass [TBI [191 EH EI] . 

Benchmark point BM1 was chosen so that A\ decays mostly to r's due to the phase space 
suppression of A\ decays to bottom quarks. In the R-axion limit m A o oc y/s, so we raised the 
mass of A\ by increasing s which leads to BM2. In BM2, A\ the branching ratio to bottom 
quarks is 0.9 while that to r-leptons is 0.06 because m A o ^> 2mb. Finally, the benchmark 
point BM3 was chosen so to illustrate that regions of NMSSM parameter space exist where 
A\ need not be light and the branching ratio of — > ZA\ can still be enhanced. These 
three specific benchmark scenarios serve as an important existence proof that it is possible 
to have significantly enhanced decay rates of — > ZA\ compared to the — > bb and still 
have a H® state with similar branching ratios as a Standard Model Higgs. However, it is 
important to keep in mind that even more optimistic scenarios are possible for the signal 
rate. For example, the A decays into taus could be enhanced in all parts of parameter space 
in a leptophilic two Higgs doublet model. 
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Figure 1: Our r reconstruction efficiency for a sample of Drell-Yan Z — > tt events where 
one tau decays in a one-prong hadronic decay as a function of the tau's generated px- The 
efficiency asymptotes to ~ 60% at high p T and we have plotted an analytic function to guide 
the eye. 

3 Limits on the simplified model from current searches 

In this section we take the benchmark points shown in Sec. [2]to be indicative of a generalized 
model where a nonstandard Higgs boson H dominantly decays into the Z boson and a 
lighter pseudo-scalar A. From now on, we proceed model-independently and analyze the 
phenomenology of the signal process H — > ZA — > Zt + t~ in the multilepton decay channel 
for a broad range of H, A masses. We start by finding the efficiency of observing this model 
in the CMS multilepton analyses under the selection cuts in Ref. [7]. 



Event simulation 

We generate samples of signal events for a broad range of H, A masses using Pythia8.170 [22] 
including the effects of initial state radiation, final state radiation, multiple interactions and 
fragmentation. These events samples were generated for pp collisions at a/s = 7 TeV using 
the CTEQ5L parton distribution functions [23] events we only consider the leptonic decay 
of the Z-boson (including r's) and assume that the A decays only into r-leptons. We do not 
apply any detector simulation on the events, in order to apply our own tau reconstruction. 

The selection criteria for the CMS multilepton analysis of Ref. [7J are as follows. The 
two trigger requirements were — single lepton: one muon (electron) has a Pt > 35 (85) GeV 
and dilepton: two of the leptons have p\ > 20 GeV and p\ > 10 GeV. An electron or 
muon is identified, if its px > 8 GeV and \r]\ < 2.1. Furthermore an isolation criterion 
/Rei = E cone /Ei < 0.15 is imposed on the leptons, where Ei is the energy of the lepton and 
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Figure 2: Contours of the efficiencies in the one-prong hadronically decay r channel, 1 
one-prong hadronically decay r channel and 2 one-prong hadronically decay r channel. In 
these channels, there is an electron or muon pair with mass consistent with the Z and a 
total of four leptons (electrons, muons, and one-prong taus). The brown hatched regions in 
the upper left are plotted to avoid the kinematically squeezed region where m H — tjia < 100 
GeV. 



E COIie is the surrounding visible energy in a cone of AR = 0.3(0.4) for muons (electrons). 
For reconstructed hadronic r-leptons we consider only one-pronged hadronic decays where 
the pt of the charged hadron is required to be > 8.0 GeV. A tau is isolated if I^ e \ < 0.15 
where I^ e \ is the ratio of total energy inside an annulus with 0.1 < AR < 0.3 to the total 
energy inside a cone with AR < 0.1. To speed up the analysis, we place a lower cut on tracks 
of pt > 0.3 GeV to be considered in the analysis and have checked that tau reconstruction 
efficiencies are insensitive to this cut. 

Ref. [7] uses the CMS particle flow algorithm to identify the neutral pions and pt of the 
Th candidate, which may lead to a difference between our simulation and the CMS data. As 
a check of our tau reconstruction, we simulated one-prong tau decays in Drell-Yan Z — > tt 
and found reconstruction efficiencies as shown in Fig. [I] as a function of the generated tau 
Pt- Our asymptotic efficiency is about 60% at high pt- This can be compared with the 
published CMS tau efficiencies in figure 3 of Ref. |24j. There is a difference in presentation 
since the CMS figures are plotted with respected to the generated visible tau pt which is 
approximately 1/2 — 1/3 of the generated tau pt- However, taking this into account, the 
behavior is similar to the TANC medium tau algorithm and systematically higher than the 
efficiencies of the HPS algorithms (which peak at 50%). This gives us confidence that our 
tau reconstruction is realistic and maybe just a bit more optimistic than the true CMS 
algorithms. 

As the Z-boson can also decay to r-leptons which can further decay into e or /i, there 
is also a possibility of this process also contributing to the multilepton events in Ref. [11] 
and Ref. [7] where the invariant mass of no two leptons are within the Z mass window 
75 GeV < rri£+£- < 105 GeV. However these electrons and muons are typically soft and 
therefore such events have more difficulty with the triggering requirement. 
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In Fig. [2] we show the contours of the efficiencies for each channel, 

_ N {selected), i 

ei = N(Z (e+e-, fi+fi-, t+t-), A t+t-) ^ 

where the channels are the 4£ + Or^, 3£ + lr^ and 2£ + 2r^ channels with a reconstructed 
leptonic Z. Thus, our convention takes out the Z branching ratio into the three generations 
of leptons and the A branching ratio into r pairs, but includes the r branching ratios into the 
efficiency. The brown hatched regions in the upper left are plotted to avoid the kinematically 
squeezed region where m# — rriA < 100 GeV. In the Or^ channel, each of the r-leptons from 
the A decay have decayed leptonically via the three-body decay. Thus, the resulting leptons 
from the decay of the r's are relatively soft. This explains the efficiency improvement as 
itla increases, since the r's have larger average px and are more geometrically leading to 
better isolation. The lr^ and 2r^ efficiency curves have a similar structure to the Ort events 
because of the similar effects of boosting the A and r-leptons. However, the value of m# 
with optimal efficiency for the lr^ and 2r^ channels is different because of the different 
isolation requirements for 1-pronged hadronic r's and the greater visible px in such decays 
as compared to the leptonic case. This helps to explain the slope of the contours, as a smaller 
boost to the A is required for hadronic taus to pass the px selection. 

Using a toy model, we can get further insight into the inefficiency of the searches at low 
rriA- In our toy model, we assume that H is produced at rest, with the A particle being 
produced in the transverse direction. This A is taken to decay into two single prong r's where 
the r's decay r + — > ir + , u T , assuming p w + = p T +/2. In this case, we can very simply predict 
the visible pion kinematics as a function of the r + decay angle in the A rest frame. The px 
and AR of the charged pions are shown in Fig. [3] for the case of m#, rriA = 200, 10 GeV. 
Here, we see that the tau reconstruction requirement of a track with px > 8 GeV restricts us 
to I cos 9cm I < -6 in order to select both hadronic r's. However, as the AR figure shows, the 
isolation condition is in direct conflict, requiring | cos 9cm I ^ -6. Thus, these two conditions 
are in tension. Due to the softness of the pions, the configuration that works best for getting 
substantial px is where the two taus decay in the longitudinal direction, so that the boost 
enhances both of their transverse momenta. However, at the same time this pushes the pions 
on top of each other, worsening isolation. This tension is exacerbated with larger boosts. 
For example, as the H mass is increased, the slope of the px plots increases whereas the dip 
of the AR plot decreases. Thus, we see that the standard tau selection and reconstruction 
is inefficient for the boosted regime. For such decays, searches for boosted taus has been 
show to be efficient [25J, in particular using iV-subjettiness [26]. Thus, multilepton analyses 
should consider a boosted tau pair object as a way to recover such regions of parameter 
space. 

Signal limits using the CMS analysis 

The multilepton analysis Ref. [7] observed number of events and expected background is 
tabulated in Tab. |4} There is a slight excess, predominantly in the 2r^ channel, which 
unfortunately is not fit well by our signal. This is due to the fact that the efficiencies for the 
1t/i, 2r/j channels are similar (see Fig. |2|), which limits the amount of the excess that can be 




Figure 3: Toy model kinematics of the charged pions as a function of the decay angle of the 
t + in the rest frame of A. The pt selection cut is shaded on the left plot and the isolation 
annulus is shaded on the right plot. The masses for these plots are m#, vtia = 200, 10 GeV. 
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Table 4: The observed number of events and expected backgrounds for the three channels as 
given in Ref. [7]. In the second to last column is our derived 95% CL. limit on the number 
of signal events in each channel and the last column is the projected limit a 30 fb _1 analysis 
would have given only statistical background errors. 



explained. In order to set limits on our model, we use this data to calculate the maximum 
number of signal events at 95% C.L. in each channel. Due to the low statistics we assume 
a Poisson distribution for the number of events and therefore the probability of observing at 
least iV bs events due to signal and background is P(A^ b s |5' + B) = r(iV b s + 1, S + B) /N Q -b s \, 
where T(A, B) is the incomplete gamma function. We assume the background has a gaussian 
distribution with mean and variance a B and find the maximum allowed number of signal 
events at 95% C.L. (S'^ ax ) by solving the equation 
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(20) 



where Mb normalizes the background's gaussian distribution over the interval [0, oo). Using 
these limits on the number of signal events in each channel for each selection, we can put 
bounds on the cross-section for this process, normalized to the Standard Model Higgs pro- 
duction cross-section by defining a signal strength parameter 



Ci Max 
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(t Hsm x BTZ(Z -> /+/-) x x £ 
e, fi, t and e l is the efficiency for i th channel and £ is the integrated luminosity. 
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Figure 4: In the left figure, the contours of the 95% CL limit on the signal parameter /z 95 
is displayed, while in the right figure, the projected limit a 30 fb -1 CMS-like analysis would 
have. 



In Fig. [4] we present contours of the minimum value of /X95 for the C = 4.8 fb -1 of data. 
The limits are weaker for small values of rriA because of the lower efficiencies and they are 
stronger for small values of m# because of the larger production cross-section. Using the 
couplings in Tab. [2] and the branching ratios in Tab. [3] we see that these benchmarks have 
H ~ 0.25 for rriA = 10 GeV and fi ~ 0.023 for > 10 GeV. Therefore the benchmark 
points considered in Tab. [I] are not constrained by the present experimental data. However, 
in more optimistic scenarios where the A is leptophilic (i.e. Br(A\ — > rr) ~ 1), we see that 
the Standard Model cross section is already strongly constrained for A masses above about 
15-20 GeV. 

With improvements in statistics and in the background modeling, a future multilepton 
analysis would have improved sensitivities to this signal and start to constrain more inter- 
esting signal rates. To estimate this improvement, in Tab. |4| we have projected the signal 
events allowed in a 30 ftr 1 CMS-like analysis. In making this estimate, we have ignored the 
background and signal cross section changes with = 8 TeV running and further assumed 
that the background uncertainty can be reduced to being purely statistical. As shown in 
the right figure of Fig. |4j this gives a projected limit at 30 fb -1 that is roughly 4 — 5 times 
stronger than the current analysis. This is very close to being sensitive to our benchmark 
signal rates and would place stringent constraints on scenarios where the A branching ratio 
to taus is enhanced. Hence, a multilepton analysis using the full LHC data set of 2012 could 
have an interesting reach for for this nonstandard Higgs signal. 
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4 Comparison of methods of mass reconstruction H — > 

ZA 



If an excess in these channels is seen in future multilepton analyses, it will be important 
to reconstruct the signal in order to determine the underlying theory. As a step in this 
direction, in this section we consider reconstructing this signal by measuring the masses m# 
and rriA through a variety of techniques. In particular we consider three possibilities: i) 
transverse mass, ii) collinear mass and Hi) an analytic solution based on trial masses for 
H, A. In this section we will be concentrating purely on the 2r/j channel because there are 
more neutrino final states in the Or^ and lr^ channels, complicating the reconstruction. We 
summarize these mass reconstruction methods below. 

Transverse Mass Variables: Using the visible components of the r's pv 12 -> the recon- 
structed Z momentum pz and the total MET components we can define the transverse 
masses 



m T A = ^p V + 2(E v El-p v .pT) (22) 
m T H = \/(pv + pzY + 2((E V + E z )El - {p T v + pj) ■ pj) (23) 



where pv = Pv^ + Pv 2 - These variables have the property that m x < mx, so measuring the 
endpoints gives a determination of the masses. 

Collinear Approximation: Since the pseudo-scalar A typically has a large boost, it is a 
good approximation to assume that the final state neutrinos are collinear with the visible 
final state hadrons of the taus. In this approximation, the neutrino momenta are proportional 
to the visible components of the r-leptons and we need to solve the linear equation [28] 

^iPv, + ^Pv 2 =Pl- (24) 

Physical solutions require the coefficients Ai,2 to be positive. From these, one approximates 
the r momenta as (1 + \i)PVi to determine the A, H masses. 

Analytic Solution: Finally, as shown in Appendix [A] for trial masses mu and itla we 
can solve the the neutrino momenta exactly. Similar to Ref. [29], for each event, there is 
an allowed region of masses where there are consistent neutrino solutions. For each such 
event, our estimator for the masses is the center of mass of the allowed (m#, m^) region. 
See Appendix |Aj for more details on this method. 

In Fig. [5]we show the histograms for each of these different mass variables for signal events 
generated with (m H , m A ) = (200 GeV, 10 GeV), (200 GeV, 100 GeV) and (300 GeV, 100 GeV). 
The blue (darker grey) histograms correspond to method iii), the average values of m# 
and itia that provide analytic neutrino solutions for each event. The yellow (lighter grey) 



histograms correspond to the transverse masses for each event as defined in Eq. (22) and 



Eq. (23). The purple (grey) histograms are the m# and rriA masses calculated in the collinear 
approximation. The peaks of the blue and purple histograms determine the reconstructed 
mass of nin and tua, while the transverse masses (shown in yellow) give another estimate of 
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Figure 5: Comparison of the reconstruction of mass variables for a (m^m^) = 
(200 GeV, 10 GeV), (200 GeV, 100 GeV) and (300 GeV, 100 GeV). The left panels are 
the reconstruction of the H mass while the right panels correspond to those of the A mass. 
The yellow (light grey) histogram corresponds to the transverse masses of H and A method 
i) and its drop off determines the mass, the purple (grey) histograms correspond to the 
mass reconstructed by using the collinear solution to the neutrino momenta method ii) and 
the blue (dark grey) histograms corresponds to the reconstructed mass using the analytic 
method in). 
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the masses as the transverse mass should drop above the physical mass. Hence the peaks for 
blue and purple histograms are to be compared with the value of the mass where the yellow 
histograms drops down. As can be seen by the plots, this drop off is not always sharp, so 
an accurate extraction of the mass from the transverse mass distribution can require a good 
understanding of its shape. 

The average solutions shown by the blue histogram typically reconstruct an m# that 
agrees well with the true value in the signal, although the value of ttla is slightly higher than 
its true value for low values of m^- This slight bias of reconstructed A masses, at low values 
of rriA, is due to the selection requirements preferring taus with a higher p?. Furthermore 
about 90% of the 2rh events that pass the isolation and trigger cuts, reconstruct to a physical 
solution. This is to be compared to the collinear solution where less than 50% of the 2th 
events lead to physical solutions with positive Ai^- The plotted event counts reflects this 
difference in solution efficiency, giving a visual indication of the increase in reconstructible 
events. Thus, for a given set of events, the analytic method can provide a more accurate and 
efficient measurement of the H and A masses. The reconstructed masses using the collinear 
method are typically slightly lower than physical masses because that angles between the 
neutrinos and visible hadrons are neglected. From Fig. |5j we see that the widths of the 
distributions using methods ii) and Hi) are comparable. Comparing the transverse mass 
distributions for itlh = 200 GeV and = 300 GeV, we see that the mass resolution of H 
and A using the end points of the distributions is better for lower masses than for higher 
masses. Therefore, the mass determination using the transverse mass distribution is more 
complicated for these higher mass points. 

To summarize, we have compared three different mass variables as a means to reconstruct 
this signal process. These mass reconstruction methods all work reasonably well, displaying 
different advantages and disadvantages depending on the situation. We focused on the 2r h 
channel, since there the kinematics can be solved once trial masses for A, H are given. 
Even though the analytic method can be quite sensitive to the measured visible momenta, 
most of the 2r^ solutions lead to a physical solution. This result is to be compared to the 
collinear method where a majority of the events do not have have a physical solution and 
the transverse mass distribution whose endpoint may be uncertain. In our reconstruction, 
the only uncertainty was the decay widths of the heavy particles, thus our analysis gives an 
idea of the irreducible uncertainties without considering detector effects. To improve further, 
there are more advanced techniques being developed which could yield further improvements. 
Likelihood methods taking into account the r decay kinematics, such as [30], would be worth 
applying to this process. 



5 Conclusion 

In this article we have shown that it is possible for nonstandard CP-even Higgs bosons to 
have large branching fraction into ZA where ^ 100 GeV and A has a sizeable branching 
fraction into r-leptons. In particular, in the NMSSM we have presented three benchmark 
scenarios where the — > ZA\ decay mode has the largest branching fraction while H\ 
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phenomenology is SM-like. In these scenarios, low tan/3 and large mixing in the pseudo- 
scalar sector are preferred. These rates can be further enhanced in more optimistic scenarios 
like in leptophilic two Higgs doublet models. 

In the collider study of the H — y ZA — y Zt + t~ scenario we found the efficiencies for 
passing the selection cuts in the Or^, Ith and 2th channels that included a leptonic Z in the 
latest CMS multilepton analysis [7j. The shape of the efficiencies contours for each channel 
are due an interplay between the selection cuts, the isolation requirements and the kinematics 
of the events. Using a toy model, we demonstrated that the low parameter space is 
inefficiently picked up by the multilepton selection. In particular, the transverse momenta 
and isolation requirements for the taus are in direct tension with each other. This motivates 
including boosted di-tau jets, as explored in [22], as a physics object in future multilepton 
analyses. Using the observed and expected background events we constructed the 95% C.L. 
limits on the number of signal events in each channel. We found that the strongest limit 
on such a scenario was due to the lr^ channel because of a high signal efficiency and the 
consistency between expected and observed events. In addition, we made a projection of 
the reach of a 30 fb _1 CMS-like study would have on this scenario. If the background 
uncertainties can be reduced to purely statistical, we estimate that a large portion of the 
interesting parameter space can be covered. 

Finally we considered the possibility of observing this scenario by measuring the m# and 
rriA masses in the transverse mass, the collinear mass and an analytical solution based on the 
trial masses in the 2rh channel. We found that this analytical solution allows for a greater 
number of events with physical solutions than the collinear method, while maintaining a 
similar resolution. The analytical method can also reconstruct larger values of the H mass, 
which may be problematic using the transverse mass distribution, as the fall off of the 
transverse mass distribution is much softer for heavier masses. 

Note Added: A new CMS multilepton analysis was recently presented at the HCP 2012 
conference using 9.2 fb~ l of 8 TeV data f31^ . The changes in analysis have not been published, 
but due to the reduction in expected number of events, the cuts have been increased. If px 
requirements have been increased, this will have a strong impact on the efficiency for our 
signal. 

A Neutrino Solution Using Trial Values of m# and tra 

In this appendix, we show how to solve for the two neutrino four vectors p ui 2 , knowing the 
momenta of the visible decay products of the r's pv 12 an d the reconstructed Z vector p z . 
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There are eight kinematic constraints 

A = = Pl 2 (25) 

(Pi* + PvJ 2 = m r = (P»a + Pv 2 ) 2 (26) 

m A = (Pvi + VVx + Pu2 + Pv 2 f (27) 

m n = (Pz + Vvx + Pvi + Pu 2 + Pv 2 f (28) 

(29) 

< + < = V\ (30) 

where are the trial masses and for simplicity we have defined 

p±=p Ul ±p U2 . (31) 

In particular, the x, y components of p + are the missing transverse energy components. 
The undetermined components of p + satisfy 

P% = EzE+ ~ Ah (32) 



Pz 



2 / E Z E + - A H \ 2 f E Z E + - A H 



El - [^r^ ) + * [E + E V - -^r^Pv) = A, (33) 
where 

Pv = PVi + Pv 2 (34) 

A H = - (m 2 H -m\- m 2 A ) - p z ■ Pv + Pz • P+ (35) 

A A = m 2 A -p 2 v + {p T + f + 2p T + -pl. (36) 

Therefore for a particular (m^,m^), a physical solution is possible only if a real positive 
root to Eq. (33) exists. 

The p- equations are 

P--p Vl =A 1 (37) 

p^-p V2 = -A 2 (38) 

P--p + = (39) 

P- = -Pi (40) 

where the terms on the right hand side are 

A; = m 2 - p\. - p + ■ p Vi (41) 

We can solve for the spatial components in terms of E_ 

p_ = E-A + B (42) 
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where 



.4 



B 



X 




Pv! Pvi Vv x 

Pv 2 Pv 2 Pv 2 
p\ p\ p\ 



(43) 
(44) 
(45) 



Plugging this solution back into Eq. (40) we see that a solution for the trial m# and 
values exists only if E_ is real. Hence the conditions for a particular event arising from a 



trial rriHifnA are that at least one of the roots of Eq. (33) is positive and real and the p + 
corresponding to that root also satisfies the condition 

(Ap + )-B(p + )) 2 -(l-\A(p + )f)(pi-\B( P+ )f) > 



which ensures that E_ is real. 
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